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We study theoretically diffusion of one-dimensional Frenkel excitons in J-aggregates at temper- 
atures that are smaller or of the order of the J-band width. We consider an aggregate as an open 
linear chain with uncorrelated on-site (diagonal) disorder that localizes the exciton at chain seg- 
ments of size smaller than the full chain length. The exciton diffusion over the localization segments 
is considered as incoherent hopping. The diffusion is probed by the exciton fluorescence quenching 
which is due to the presence of point traps in the aggregate. The rate equation for populations of 
the localized exciton states is used to describe the exciton diffusion and trapping. We show that 
there exist two regimes of the exciton diffusion at low temperatures. The first, slower one, involves 
only the states of the very tail of the density of states, while the second, much faster one, also 
involves the higher states that are close to the bottom of the exciton band. The activation energy 
for the first regime of diffusion is of the order of one fifth of the J-band width, while for the second 
one it is of the order of the full J-band width. We discuss also the experimental data on the fast 
low-temperature exciton-exciton annihilation reported recently by I. G. Scheblykin et al, J. Phys. 
Chem. B 104, 10949 (2000). 
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I. INTRODUCTION 



Since the seminal works by JelleyQ I and Scheibe Q , 
the concept of Frenkel excitons H 0> 3 nas been used 
for explaining of the remarkable optical properties of 
molecular J-aggregates: (i) the appearance of a nar- 
row and intense line in the red-wing of the absorption 
spectra (so called J-band), the full width of which is 
of the order of several tens of wavenumbers at cryo- 
genic temperatures and (ii) the increase of the oscil- 
lator strength of the optical transition by almost two 
orders of magnitude 0, 0, IS 01 fioj - During the 
nineties, a considerable progress in understanding of lin- 
ear and nonlinear optical dynamics of J-aggreg ates was 
made (for details see the reviews Il2t Il3l and ref- 
erences therein). In spite of the fact that monomers 
which form the aggregates have complex chemical struc- 
ture, both linear and nonlinear optical dynamics in J- 
aggregates have been successfully described on the basis 
of the simplest one-dimensional (ID) tight-binding model 
with diagonal and/or off-diagonal disorder, both uncor- 
related pi [ii El and correlated El E [H 111 ■ 

The eigenstates of a homogeneous (non-disordered) 
J-aggregate extend over the whole (N monomers) ag- 
gregate. Disorder localizes the lowest in energy exci- 
ton states at segments of about N* molecules; N* de- 
pends on the disorder magnitude and is typically much 
smaller than the total number of molecules in the chain: 
N* -C N. One of the most important consequences of 
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this localization is the appearance of states below the 
bottom of the bare exciton band. These states form the 
tail of the density of states (DOS) and carry almost the 
whole oscillator strength of the aggregate. For this rea- 
son the one-exciton absorption in J-aggregates is spec- 
trally located at the tail of the DOS (see, for instance, 
Refs. 0, llOri and the width of the absorption band is of 
the order of the width of the DOS tail. 

The exciton diffusion in a disordered aggregate is es- 
sentially the transition from one localized eigenstate to 
another. The transition probability depends particularly 
on the temperature, the energy spacing between the in- 
volved states and the overlap of these states. The lower 
states, being localized at different A r *-molecule segments 
of the aggregate, overlap very weakly Contrary to 
that, the higher exciton states, that are localized at seg- 
ments larger than N* molecules, overlap strongly with 
several lower tail states. Although higher states are ther- 
mally less favorable, the hops from the lower to higher 
states can be faster than between the lower states because 
of the higher overlap. In this paper, we show that this 
competition between the overlaps and thermal favorabil- 
ity result in a complex scenario of the exciton transport 
at low temperatures. At zero temperature an exciton re- 
sides in one of the lower states at the tail of the DOS. 
As temperature rises, first, the exciton starts to diffuse 
over the weakly overlapped states of the DOS tail. The 
activation energy for this regime is of the order of 1 /4 of 
the DOS tail width (that is of the order of the J-band 
width). The diffusion in this regime is very slow. As 
the temperature increases further, the higher states come 
into play. As these states overlap much better with the 
lower states and each other and also are more extended, 
the diffusion rate increases by several orders of magni- 
tude. The activation energy for this faster regime of the 
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exciton diffusion is of the order of the DOS tail width or, 
in other words, of the order the J-band width. 

To the best of our knowledge, these aspects of the ID 
diffusion problem have not been discussed in the litera- 
ture yet. The same tight-binding Hamiltonian was used 
to describe transport properties of electrons in doped 
semiconductors [2(| as well as those of optical excita- 
tions in activated glasses 0, |U . It is to be stressed that 
despite the seeming similarity of these problems the out- 
lined scenario of the low temperature ID diffusion over 
the localization segments is more complex than the dif- 
fusion over one-level point impurity centers in semicon- 
ductors or glasses. The major complication comes from 
the fact that the exciton can hop sideways to a differ- 
ent segment not only directly (a "horizontal" hop) but 
also indirectly via higher states, so that "vertical" hops 
up in energy become extremely important. In fact, it is 
the indirect hops that provide the dominant contribution 
to the diffusion rate at temperatures of the order of the 
J-band width. 

We use the quenching of the exciton fluorescence by 
point traps to probe the exciton diffusion. The temper- 
ature range lower or of the order of magnitude of the 
J-band width is of our primary interest; higher tempera- 
tures are beyond the scope of the present work. 

The outline of the paper is as follows. In Sec. [H] we 
present the microscopic model of exciton trapping. Sec- 
tion IIIII is focused on the qualitative discussion of the 
channels of the exciton diffusion over the localization 
segments. The results of our numerical simulations of 
the exciton fluorescence quenching, obtained on the ba- 
sis of the rate equation approach, are the contents of 
Sec. II VI In Sec. \V\ we conclude the paper and discuss 
the results of the recent experiments on the fast exciton- 
exciton annihilation in the aggregates of the triethylth- 
iacarbocyanine salt of 3,3'-bis(sulfopropyl)-5,5'-dichloro- 
9-ethylthiacarbocyanine (THIATS) 

II. MICROSCOPIC MODEL OF THE EXCITON 
FLUORESCENCE QUENCHING 

We model a J-aggregate by N (N 3> 1) optically ac- 
tive two-level molecules forming a regular in space ID 
open chain. The corresponding Frenkel exciton Hamilto- 
nian reads Q (for the sake of simplicity only the nearest- 
neighbor interaction is considered) 

N JV-1 

H = Y / E n \n)(n\-jY / (\n+l)(n\ + \n)(n+l\) . (1) 

n— 1 n— 1 

Here E n is the excitation energy of the n-th molecule, 
\n) denotes the state vector of the n-th excited molecule. 
The energies E n are assumed to be Gaussian uncorre- 
cted (for different sites) stochastic variables distributed 
around the mean value u>o (which is set to zero without 
loss of generality) with the standard deviation A. The 
hopping integral, — J, is considered to be non-random 



and negative ( J > 0), which corresponds to the case of J- 
aggregates (see, e.g., Ref. @). In this case the states cou- 
pled to the light are those close to the bottom of the ex- 
citon band. In what follows, moderate disorder (A < J) 
is considered. This implies that the exciton eigenstates 
ip v (v = 1, 2, . . . , N), found from the Schroedinger equa- 
tion 

N 

H nm (p vm = e v (p vn , H nm = (n\H\m) , (2) 

m— 1 

are extended over relatively large segments of the chain. 
However, the typical size of these localization segments, 
N* , is small compared to the full chain length N (units 
of the lattice constant are used throughout the paper). 

Having been excited into an eigenstate is, an exciton 
cannot hop to other eigenstates if coupling to vibrations 
is not taken into account. We assume that this coupling 
is weak and do not consider polaron effects. This limit 
is applicable to a number of J-aggregates as the Stokes 
shift of the luminescence spectra with respect to the ab- 
sorption spectra is usually small 0, 0- The exciton- 
vibration interaction causes the incoherent hopping of 
excitons from one eigenstate to another. We take the 
hopping rate from the state v to the state \i in the fol- 
lowing form (see, e.g., Ref. 0) 

JV 
n=l 

1 1 + n(e u - Efj), < e v 

Here, the constant Wq characterizes the amplitude of 
hopping and n(e) = [exp(e/T) — I] -1 is the occupa- 
tion number of the vibration mode with the energy 
e (the Boltzmann constant is set to unity). Due to 
the presence of the n{e) and I + n(e) factors, the rate 
W^ v meets the principle of detailed balance: — 
W V fj,exp\{e v — e^)/T]. Thus, in the absence of decay 
channels, the eventual exciton distribution is the Boltz- 
mann equilibrium distribution. The sum over sites in @ 
represents the overlap integral of exciton probabilities for 
the states fi and v. The spectral factor S{\e v — e^l) de- 
pends on the details of the exciton-phonon coupling as 
well as on the DOS of the medium into which the ag- 
gregate is embedded. For example, within the Debye 
model for the density of phonon states, this factor takes 
the form S{E V - E») = {\E V - E^/Jf 25]. However, 
this model is applicable to glassy media (the media we 
assume as the host) only in a narrow frequency interval 
of the order of several wavenumbers (see, for instance, 
Ref. |M Therefore, as in Refs. HI HI, we re- 

strict ourselves to a linear approximation to this factor, 
S(E V — Efj) = \E V — Epl/J. This accounts for reduc- 
tion of the exciton- vibration interaction in the long-wave 
acoustic limit 0,0- Also, it eliminates the divergence of 
W ViL at small values of \E V — E^\. 
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The diffusion of Frenkel excitons can be probed by 
quenching of the exciton fluorescence by traps. Con- 
sider an aggregate with point traps, namely monomers 
at which an exciton decays non-radiatively and very fast 
compared to the typical spontaneous emission rate of the 
aggregate. Then those excitons that reach the traps de- 
cay non-radiatively and contribute to the fluorescence 
quenching. If an exciton is created far from the trap 
it has to diffuse to the trap to be quenched, the faster it 
diffuses the more effective is the fluorescence quenching. 
Thus, the quenching rate depends on the diffusion rate 
and can be used as a probe of the latter. 

We define the quenching rate of the exciton state v as: 



where 



N„ 



(4) 



where T is the amplitude of exciton quenching and the 
sum runs over positions of the N q traps. Thus, we take 
the quenching rate to be proportional to the probability 
to find the exciton at trap sites. 

We describe the process of the exciton trapping by 
means of the rate equation: 



Pu 



N 

E 



{lu + T V )P V + V(WV P„ - P„) 



(5) 



where P v is the population of the z/th exciton eigenstate 
and the dot denotes the time derivative, 7^ =7/1/ is the 
spontaneous emission rate of the i/th exciton state, while 
7 is that of a monomer, /„ = (X) n =i ¥W) 2 being the 
oscillator strength of the state v. 

The temperature dependence of the exciton quenching 
is calculated as follows. We admit the definition of the 
exciton fluorescence decay time as the integrated total 
population p5|: 



dt 



N 

£■ 

i/=i 



(0) 



where angle brackets denote the average over disorder 
realizations and traps positions. The decay time has to 
be calculated for aggregates with traps (denoted as r) 
and without traps (denoted as To). The quenching rate 
is then defined as 



N„=0 



1 

7-Q 



(7) 



This quantity carries information about the diffusion rate 
and is the object of our analysis. 

The definition of the decay rate as the integrated total 
population allows for considerable simplification of the 
calculation procedure. Write the solution of Eq. J5J in 
the formal matrix form 



N 

PAt) = E ( £ 



-Ri 



H=l 



PM 



(8) 
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+ E W ^ <V - W vlt . (9) 



After substitution of © into Eq. © and integration over 
time, t can be expressed in terms of the i?-matrix: 



JV 



t = E (^ 1p m(o) 

1 



(10) 



Calculation of the quenching rate W q requires the cal- 
culation of the inverse matrix for each realiza- 
tion of disorder rather than the fluorescence kinetics. 
The inverse matrix is to be found twice: for an ag- 
gregate with and without traps. Note that the decay 
time Tn also depends on temperature (see, for example, 

Refs. HIES H|). 



III. QUALITATIVE ARGUMENTS 

At low temperatures excitons reside in the tail of the 
DOS, that is, below the bottom of the bare exciton band, 
E = — 2 J. As we show below, higher states that are 
close to the bottom of the bare band contribute to the 
exciton diffusion as well. Therefore, these two parts of 
the exciton energy spectrum are of primary importance 
for the low-temperature exciton transport. 



A. Analyzing the low energy structure 

Here, we recall briefly the concept of the local (hid- 
den) energy structure of localized ID excitons [rRl3lIl32T |. 
which was proved to exist in the vicinity of the band bot- 
tom [lJiyJ]- According to this concept, the low-energy 
one-exciton eigenfunctions obtained for a fixed realiza- 
tion of disorder are localized at segments of typical size 
N* (localization length). Some of these localized states 
(about 30%) can be grouped into local manifolds of two 
(or sometimes more) states that are localized at the same 
A r *-molecule segment (see the states filled with black 
color and joined by ellipses in Fig.^l. It turns out that 
the structure of the exciton states in each local manifold 
is very similar to the structure of the lower states of a 
homogeneous (non-disordered) linear chain of length N* . 
In particular, the lowest state in a manifold has a wave 
function without nodes within its localization segment. 
Such a state can be interpreted as the local ground state 
of the segment (italic is used to distinguish this state 
from the true ground state, that is, the state with the 
lowest energy in each realization). A local ground state 
carries large oscillator strength, approximately N* times 
larger than that of a monomer, so that the typical spon- 
taneous emission rate is 7* = 7 N* . The scaling law of 
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FIG. 1: The energy structure of the exciton states in the 
vicinity of the bottom of the exciton band. The states are 
obtained by diagonalization of the Hamiltonian Q for a linear 
chain of 300 molecules and the disorder magnitude A = 0.1J. 
The baseline of each state represents its energy in units of the 
spacing in the local energy structure £12. The origin of the 
exciton energy is set to the lowest energy for the realization. 
The wave functions are in arbitrary units. It is clearly seen 
that some lower states can be grouped into local manifolds 
(the manifolds are joined by ellipses) . The states within each 
manifold are localized at the same segment of typical length 
N* (this length is given by the bar in the lower right corner) , 
they overlap well with each other and overlap much weaker 
with the states localized at other segments. The higher states 
(filled with gray color) are more extended than the lower ones. 
Typically, they overlap well with several lower states and with 
each other. 



the localization length is |3j 



A = S.7 I j 



-0.67 



(11) 



The energy distribution of the local ground states, cal- 
culated as described in Ref. |3J], and the total DOS are 
presented in FigEl This figure shows that almost all local 
ground states belong to the tail of the DOS, as has been 
mentioned in the Introduction. 

The second state in a manifold has a node within the 
localization segment (see the states filled with black color 
and joined by ellipses in Fig.^l and looks like the first lo- 
cal excited state of the segment. Its oscillator strength is 
typically an order of magnitude smaller than that of the 
local ground state. It is important to recall here that, con- 
trary to the eigenstates from the same manifold, the lower 
states localized at different segments overlap weakly (see 
all states filled with black color in Fig.^). The energies 
of the local ground states are distributed within the inter- 
val \f2o-\i ((T11 being the average spacing between local 
ground states). This interval is larger than the typical 
ener gy spacing £12 between the levels in a local mani- 
fold J34J. For this reason, the local energy structure can- 
not be seen either in the DOS (see Fig. or m the linear 
absorption spectra (see, for instance, Ref. 10]). However, 
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FIG. 2: The total DOS (solid line) and the DOS of the local 
ground states P{E\) (dashed line) for different magnitudes 
of disorder A. The DOS is normalized to N, P(Ei) is nor- 
malized to N/N* . The vertical lines indicate positions of the 
curves' maxima. For all considered magnitudes of disorder 
the maximum of the local DOS is shifted with respect to the 
maximum of the total DOS by about the mean spacing in the 
local energy structure, £12. 



this structure determines the nonlinear optical response 
of the aggregate |3^, |3|| . 

Higher states are more extended than the local states 
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as the localization length increases with energy (see the 
states filled with gray color in Fig. Therefore, the 
higher states cannot be included into any particular local 
manifold: their wave functions covers more than one N*- 
molecule segments. Nevertheless, as all these states are 
close to the maximum of the DOS (see Fig.|2J), the typical 
energy spacing between the higher states and the covered 
local states is of the order of £12. Thus, the energy £12 
is expected to be the characteristic energy of the exciton 
diffusion. 

It is clear from the above arguments that at temper- 
atures T < e 12 , it is the states from the local mani- 
folds that determine the exciton diffusion. Two types 
of hopping over these states can be distinguished: intra- 
segment hopping and inter-segment one, involving the 
states of the same local manifold and those of different 
manifolds, respectively. As the states from different local 
manifolds overlap weakly, only inter-segment hops to ad- 
jacent segments are of importance. The disorder scaling 
of the overlap integrals — J2 n <ffc n <f% n for the local 
states of the same and adjacent segments was obtained 
in RefiEi 

/AX ' 70 

hi = 0.14 - (12a) 



B. Hopping at zero temperature 

At zero temperature an exciton can hop only down to 
lower states. Let us assume that it is in the local excited 
state 2. Then it can either hop to the local ground state 1 
of the same segment or to a lower state v' localized at an 
adjacent segment (see Fig. |3 T = 0). Because the intra- 
segment hopping is faster than the inter-segment one, 
first, the exciton hops down to the local ground state 
1 with the typical energy loss £12 (£12 being the mean 
energy spacing in the local energy structure, see Fig. |3J 
T = 0). From the local ground state, the exciton can 
hop only to a state v 1 of an adjacent segment, provided 




/ A X0.75 

I*\ » Iw2 = 0.0025 - . (12b) 



Hereafter, the indices 1 and 2 label the local states of 
the same segment while those with primes label the local 
states of an adjacent segment. As follows from Eq. 1)12(1. 
the intra-segment overlap integral is typically two orders 
of magnitude larger than the inter-segment one. Note, 
that both overlap integrals scale approximately propor- 
tional to the inverted N* (compare Eqs. I|12a|) - (|12b|l with 
Eq. (|llfl ). This proportionality holds for two exponential 
functions extended over the length N* and separated by 
the distance of the same order of magnitude, N*. 

The intra-segment hops do not result in the spatial 
displacement of excitons. Only the inter-segment hop- 
ping gives rise to the spatial motion. Nevertheless, we 
show below that both types of hops are important for 
understanding the features of the low-temperature exci- 
ton transport. 

The overlap integrals between the local states of a seg- 
ment and the higher states which are extended over this 
segment and a few adjacent ones (see the states filled 
with gray color in Fig. ^) are of the order of Ii2- This 
fact implies that even at T < £12, the indirect hops via 
these higher states can be more efficient than the direct 
inter-segment hops over the states of local manifolds (see 
below). Our calculations support this assumption. 



FIG. 3: Schematic view of exciton hopping at zero and non- 
zero temperatures. The indices 1 and 2 label the local ground 
and the first local excited states of the same segment, respec- 
tively. The v' state is localized at an adjacent segment. The 
index 3 label a higher state, which extends over two adja- 
cent segments. Hops are shown by straight arrows; the arrow 
thickness represents magnitude of the corespondent hopping 
rate. Thin wavy arrows show spontaneous emission. Its rate 
is the slowest among all, which corresponds to the limit of 
fast diffusion. 

£ v ' < £1 and the spontaneous emission rate of the local 
ground state 71 is small compared to the intra-segment 
hopping rate H^'iIt^- Hereafter, such a relationship 
between these rates is referred to as the limit of fast dif- 
fusion; only this limit is considered in this work. The 
typical energy loss during such sideways hop is of the or- 
der of the average spacing between local ground states, 
o~n (en is of the order of the J-band width). Thus, al- 
ready after one such sideways hop the exciton resides in 
a state in the tail of the DOS (see Fig. [2J. Therefore, the 
number of states with even lower energies decreases dra- 
matically, which results in a strong increase of the typical 
distance to those states and decrease of the probability to 
hop further sideways. Then the exciton either relaxes to a 
lower state of the same segment (if there is one) or decays 
spontaneously, i.e. this type of the spatio-energetic dif- 
fusion (towards lower energies) stops very quickly. Note 
that this diffusion would manifest itself in the red shift of 
the exciton emission spectrum relative to the absorption 
spectrum. The experimental data shows that such red 
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shift is either absent [g, l?J or is smaller than the J-band 
width [23L l3S| . These experimental facts indicate unam- 
biguously that at low temperatures, T <§C a%i, excitons 
make few hops before they decay due to the spontaneous 
emission, as was argued in Refs.|3Jj|40|,|4l|. Consequently, 
the zero-temperature exciton quenching is expected to be 
weak provided the concentration of quenchers is low, the 
case we are interested in. 



C. Hopping at non-zero temperatures 

At non-zero but low temperatures (0 < T < £12), an 
exciton can also hop up in energy. Consider an exciton 
in one of the lower states in the tail of the DOS, e.g. in 
the local ground state 1 (see Fig. | T > 0). For the 
reasons discussed above, first, the exciton hops up to the 
first local excited state 2 of the same segment, provided 
the hopping rate for the considered temperature is larger 
than the spontaneous emission rate 71 of the initial state 
1. During this process the exciton typically gains the 
energy £i 2 . As £12 is of the order of an j^J], already 
after the first hop up the exciton leaves the tail of the 
DOS (see Fig. |5J) and, hence, it is likely to have a lower 
state v' localized at an adjacent segment. A hop down 
to this state with loss in energy is favorable and results 
in the spatial displacement of the exciton, i.e., in the 
exciton diffusion. We stress that although only sideways 
hops result in the spatial displacement of the exciton, it 
is the initial hop up from the local ground state 1 to the 
local excited state 2 that triggers the diffusion. 

Another way for the exciton to hop sideways to the 
state v' is via the higher state 3 that overlaps well with 
both states 2 and v 1 (see Fig. [31 T > 0). As it has been 
mentioned, such hops compete with the sideways hops 
over the local states; although the hop up to the state 
3 is thermally unfavorably, the overlap integral for this 
hop, Z31, is large compared to that for an inter-segment 
hop, I y i\ . We show later that this channel of the diffusion 
becomes dominant even at relatively low temperature. 



IV. TEMPERATURE DEPENDENCE OF THE 
QUENCHING RATE 

In this section, we discuss the results of numerical cal- 
culation of the quenching rate W q . We consider the initial 
condition where the leftmost local ground state is excited 
while a single trap is located in the center of the localiza- 
tion segment of the rightmost local ground state. In this 
case, the exciton quenching is most affected by the diffu- 
sion, as the created exciton has to travel over almost the 
whole chain to be quenched. Thus, the exciton quench- 
ing at low concentration of traps can be studied. The 
quenching rate was calculated as described in section ITT1 
for the parameter set corresponding to the limit of fast 
diffusion and fast quenching (the latter limit is defined 
below) . 



A. Numerical results 

As it was already mentioned in Sec. II I II in the limit 
of fast diffusion the inter-segment down-hopping rate is 
large compared to the typical spontaneous emission rate 
of a local ground state: 

W„/i| T =o ~ W ^j- hn » 7 r. (13) 

If a quencher is located within the localization segment 
of a local state then the typical quenching rate for this 
state is T* — T/N* [see Eq. Q]. As we are interested 
in the limit of fast quenching, this rate should be taken 
larger than the typical intra-segment down-hopping rate: 

£1 9 

T/N* » W 12 \ T =o ~ Wo -j- hi . (14) 

This ensures that once an exciton hops to a local state 
of the segment with the trap, it is quenched almost in- 
stantly. 

The inequalities (|13|) and H14[) yield the relationship 
between the rate equation parameters in the limit of fast 
diffusion and quenching: 

^o^^»7 d5a) 



Wo -j In N* < T . (15b) 

The scaling laws of the values of an and £12 were ob- 
tained in Ref. [34| and are given by 

, A x 1.33 

= 0.7 J - , (16a) 



/ A x 1.36 

£12 = 0.4 J ( -J . (16b) 

On the basis of the scaling laws (|12|) and (|16fl the param- 
eter set was chosen so that for each magnitude of the dis- 
order A the following equalities hold: T* = 10VFi2|r=o 
and Wi'i|t=o = 10 7*. Calculations were performed for 
N = 1000 and averaged over 100 realizations of the dis- 
order. 

Figure 0] shows the temperature dependence of the 
quenching rate W q for different magnitudes of the disor- 
der A. In each plot, the quenching rate is given in units 
of the typical exciton radiative rate 7* = jN*. The tem- 
perature is given in units of the mean energy spacing in 
the local energy structure £12. Note that both N* and 
£12 depend on A as described by (fTT]) and (|16bl) . Figure 
^demonstrates very clearly that for all considered values 
of A at temperatures lower than £12 the quenching rate 
is vanishing. This indicates that the diffusion at these 
temperatures is not fast enough for the exciton to reach 
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Temperature, units of 8 



FIG. 4: Temperature dependence of the quenching rate W q 
calculated for a linear chain of the length N = 1000 and differ- 
ent magnitudes of the disorder A. The averaging is performed 
over 100 disorder realizations. For each realization of the dis- 
order, the leftmost local ground state is excited, while the only 
trap is located in the center of the localization segment of the 
rightmost local ground state. 



the quencher during its (spontaneous) lifetime: it emits 
a photon before it is trapped. On the contrary, just after 
the temperature exceeds approximately £12 the quench- 
ing becomes noticeable: the exciton partly diffuses to the 
trap where it decays mostly due to quenching. Specifi- 
cally, temperature of the order of 2 £12 are required for 
the quenching to become as effective as the spontaneous 
emission: W q ~ 7* = ^yN*. 



B. Discussion 

In order to understand which states contribute most 
into the quenching process it is useful to estimate the 
effective sideways hopping rate W, which is required to 
reach the quenching level W q ^7*. To do this, con- 
sider the sequence of localization segments as an effective 
chain of "sites", the typical number of which is equal to 
the number of segments, N s = N/N*; the mean spac- 
ing between these "sites" is N*. The exciton diffusion 
coefficient is then estimated as D ~ W N* 2 (the lattice 
constant is set to unity). For the quenching to be as ef- 
fective as the spontaneous decay, the exciton has to be 
at the position of the trap (located on the opposite side 
of the chain) during the lifetime 7* _1 , i.e., it has to dif- 
fuse over the distance N during this time. Equating the 
diffusion length yjD/^f* to N, we obtain the estimate for 
the required diffusion rate W: 

W-^iN/N*) 2 . (17) 

The localization length N* is equal to 38, 25 and 18 for 
A = 0.1, 0.2 and 0.3, respectively. Thus, the correspond- 
ing diffusion rates W are estimated as 6257*, I6OO7* and 
25007* (for N=1000). These values are about two orders 
of magnitude larger than the rates of sideways hops over 
the local states, taken to be IO7* in all calculations. This 
indicates that when the quenching rate becomes compa- 
rable to the spontaneous emission rate, the exciton does 
not hop directly between the local states of adjacent seg- 
ments (with the typical rate W v >\ ~ IO7*). It rather hops 
via the higher states that extend over more than one N*- 
molecule segments (see the discussion in Sec. IIII|I . The 
hopping rate via such states for T ~ 2 £12 is of the order 
of W12 which is about two orders of magnitude larger 
than W v '\ (see Eq. © and the scaling laws I|12|I % 1. 

In order to prove the above finding we performed cal- 
culations of the quenching rate W q varying the number 
of states considered in Eq. I)l()[l. More specifically, we 
considered all states up to some (variable) cut-off state. 
Figure shows the results of such study performed for 
A = 0.1J. As it can be seen from the figure, W q de- 
pends drastically on the number of states participating 
in the quenching process. In the region where W q > 7*, 
approximately 6N/N* states are required to reach the 
true value of the quenching rate that is calculated for all 
states (compare dashed and solid lines). Thus, at tem- 
peratures T > £12 the higher states provide the dominant 
contribution into the exciton quenching process. 

Figure shows the regions of the DOS that corre- 
spond to different numbers of states that were used for 
calculation of the data presented in Fig. [SJ The higher 
states lie just after the local ones, close to the maximum 
of the DOS (see Fig. EJ. Therefore, the typical energy 
spacing between the local and higher states is about £12. 
As the higher states extend over several, but not very 
many, AT*-molecule segments (see Fig.^J, the overlap in- 
tegral between these states and the covered local states is 
large. These two factors ensure high hopping rate from 
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FIG. 5: Temperature dependence of the quenching rate W q 
calculated for A = 0.1J (TV* w 40) and different numbers 
of states considered in Eq. 1101 : solid line — all TV = 1000 
states, dashed line — 6TV/TV* = 150 states, dashed-dotted 
line — 4TV/TV* = 100 states, dotted line — 3TV/TV* = 75 
states, dashed-dashed-dotted line — 2N/N* = 50 states. The 
averaging is performed over 100 disorder realizations. For 
each realization of the disorder, the leftmost local ground state 
is excited, while the only trap is located in the center of the 
localization segment of the rightmost local ground state. 



800 
600 
400 
200 



1 .3 N/N* 


i 

^ 2 N/N* 


i 

DOS 








3 N/N* 










4 N/N* 










6 N/N* 

— ■ l 













-2.1 



-2 -1.9 
Energy, units of J 



-1.7 



FIG. 6: The total DOS calculated for TV = 1000, A = 0.1J 
(TV* « 40). The DOS is normalized to TV. The vertical lines 
show the maximum energies corresponding to different num- 
bers of states considered in the rate equation that was used 
to calculate the dependencies plotted in Fig. (in the sense 
that all states lower than the specified energy are consid- 
ered). Note that the tail of the DOS is formed by 1.3TV/TV* 
states, namely, by the states of the local manifolds (N/N* 
local ground states plus 0.3TV/TV* of the local excited states; 
recall that about 30% of the local ground states form the dou- 
blets). 



the lower local to the higher states. Another important 
point is that higher states are well overlapped and more 
extended, so hops between them are typically faster and 
longer than those between the local ones. Also, the higher 



states have small oscillator strength, so as long as an exci- 
ton remains in these states it does not decay radiatively. 
The above qualitative arguments explain the dominant 
contribution of the higher states into the exciton diffusion 
and quenching within the temperature range T > E\2- 
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FIG. 7: Temperature dependence of the quenching rate W q 
calculated for A = 0.2 J (TV* » 25). Solid line — all TV = 250 
states, dashed line — 2TV/TV* = 20 states. 

It is also seen from Fig. [5] that the difference between 
the true value of the quenching rate W q and that calcu- 
lated for a restricted number of states decreases at lower 
temperatures. Figure d demonstrates the temperature 
dependence of W q obtained for A = 0.2 J (TV* ps 25) and 
temperatures T < S\2- The solid line presents the depen- 
dence calculated for all TV = 250 states considered in the 
rate equation, the dashed line — that for 2TV/TV* = 20 
states. These 2TV/TV* states include all the states of the 
local manifolds (1.3TV/TV*) and 0.7TV/TV* higher ones. 

The parameters of the rate equation was set as fol- 
lows: T* = 10Wi 2 |t=o and WVi|t=o = IOO7*. The 
chain length TV = 250 is chosen so that the effective hop- 
ping rate W ~ 7* (N/N*) 2 estimated as discussed above 
(see Eq. (|17H ) is equal to the rate of the direct hopping 
to an adjacent segment, WVi|t=o- This yields the equa- 
tion (N/N*) 2 = 100 for the chain length. This condition 
ensures that the diffusion over the lower states only can 
provide the quenching rate W q ~ 7* at T ~ £12. Indeed, 
W q fa 0.5 7* at this temperature for 2 TV/TV* states. The 
most important point demonstrated by Fig. is that be- 
low the temperature T\ pa 0.25 £12 the two curves deviate 
slightly, which means that the contribution of the higher 
states into the diffusion becomes negligible: the exciton 
hops mostly over the DOS tail states. On the contrary, 
above the temperature T\ the higher states provide the 
dominant contribution to the diffusion and quenching. 
Note also, that the value of the quenching rate at the 
critical temperature T\ is typically very small, so the ex- 
perimental observation of this "regime change" is a chal- 
lenging task. 

The critical temperature T\ at which the higher states 
come into play can be estimated by equating the typical 
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rate of the direct sideways hopping from a local state 2 to 
an adjacent local state 1' to the "vertical" hopping rate 
from the local state 2 to a higher state 3 (see Fig. T> 
0): Wi'2 — Wz2 ~ Wi2|t=o ex p(£i2/Ti)- This equation 
yields the temperature T\. 



ln(/i 2 //i' 2 ) 



£12 ~ 0.25 £i2 . 



(18) 



We stress that the numerical factor 1/ ln(ii2//i'2) ~ 0.25 
is almost independent of the disorder, as the disorder 
scalings of the overlap integrals are almost the same (see 
Eqs. l(T2"a|) - (|12bl) ). So, the estimate Ti « 0.25 e 12 is uni- 
versal for a wide range of the disorder degree. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we study theoretically the peculiarities of 
the low-temperature diffusion of the ID Frenkel excitons 
localized by a moderate diagonal disorder. The exciton 
motion over localized states is considered as incoherent 
hopping. The diffusion is probed by the exciton quench- 
ing at a point trap. We consider a single trap located 
at one end of the aggregate while the exciton is created 
initially at the other end. In this case the exciton has 
to travel over almost the whole chain to be quenched. 
Under this conditions, the quenching rate carries direct 
information about the diffusion length that the exciton 
travels over during its lifetime. The exciton quenching is 
described by the rate equation with the quenching rate 
being proportional to the probability of finding the exci- 
ton at the trap site. 

Both our qualitative arguments and numerical simu- 
lations show that there exist two regimes of the exci- 
ton diffusion. At lower temperatures, those smaller than 
T\ w 0.25 xj-band- width, the exciton diffuses mostly 
over weakly overlapped DOS tail states which determine 
the optical response and form the J-band. This regime 
of diffusion is very slow; the exciton cannot diffuse over 
large distance during its lifetime at these temperatures. 

At higher temperatures, the higher states come into 
play. The diffusion begins to built up due to the two-step 
hops via higher states. This accelerates the exciton diffu- 
sion drastically, so that an exciton can diffuse over large 
distances during its lifetime. The higher states begin to 
contribute dominantly to the diffusion at temperatures 
higher than about T\. However, the diffusion becomes 
really fast (in the sense that the quenching rate becomes 
comparable to the spontaneous emission rate of the ag- 
gregate) only at the temperatures of the order of the 
J-band width. 

In Ref. [2i| , the anomalously fast low-temperature dif- 
fusion of Frenkel excitons in J-aggregates of THIATS 
was reported. The authors of Ref. 23 studied experi- 
mentally the exciton-exciton annihilation and found that 
this effect is pronounced even at T = 5 K (3.5cm -1 ), 
while the width of J-band of THIATS J-aggregates is 
82cm -1 . In order to explain the experimental data, the 



authors assumed that an exciton travels over about 10 4 
dye molecules during its lifetime to meet another exci- 
ton and annihilate. They found also that the activation 
energy of the exciton diffusion was 15A^ (10.5cm -1 ) and 
interpreted this energy as the typical energy spacing be- 
tween the states of adjacent localization segments. 

Despite the fact that the exciton-exciton annihilation 
should to be treated differently from the exciton quench- 
ing, the model we are dealing with can easily be adapted 
for qualitative analysis of the exciton-exciton annihila- 
tion: one of the two excitons can be considered as an im- 
mobile trap for the other, while the other diffuses twice 
as fast. As reported in Ref. j2^|, the fluorescence spec- 
trum of THIATS J-aggregates is narrowed by approxi- 
mately 26cm -1 and red-shifted by 23cm -1 as compared 
to the absorption spectrum. These results indicate un- 
ambiguously that excitons make sideways hops during 
their lifetime, i.e., the rate of sideways hops over local 
states is larger than the exciton spontaneous emission 
rate. Thus, the conditions for the exciton diffusion in 
THIATS J-aggregates are similar to those studied in the 
present paper (the limit of fast diffusion). 

Discussing the above mentioned experimental data and 
its interpretation presented in Ref. ,23j, the following 
points can be made. First, the typical energy spacing 
between the states of the adjacent segments is of the or- 
der of the J-band width that is 82cm -1 and not 
10.5cm -1 . The latter value is closer to 0.25 x 82cm -1 , 
i.e. this temperature could be related to the tempera- 
ture T\ , the activation energy of the faster exciton diffu- 
sion regime. Above this temperature an exciton diffuses 
mostly over the higher states and not over the DOS tail 
states, as it was suggested in Ref. ^2^. Another point, 
and a more important one, is of the quantitative nature: 
The typical size of the localization segment in THIATS 
J-aggregates is N* — 30 42] . In the model we are using, 
this corresponds to the disorder magnitude A w 0.2 J. 
Our numerical data obtained for a chain of N = 1000 
molecules demonstrates that for this value of the disor- 
der the exciton quenching is vanishingly small for the 
temperatures T ~ (10.5/82) x £12 (we remind that £12 is 
of the order of the J-band width). In other words, the 
exciton created in the leftmost local ground state cannot 
diffuse over the whole chain of 1000 monomers during 
its lifetime. However, it can do so at the temperatures of 
the order of T ~ £12 ~ 82cm -1 . Thus, understanding the 
fast low-temperature diffusion in the aggregates of THI- 
ATS dye molecules, observed in Ref. [23, still remains 
an open question. 
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